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A maneuver design method that is particularly well-suited for determining the stability 
and control characteristics of hypersonic vehicles is described in detail. Analytical 
properties of the maneuver design are explained. The importance of these analytical 
properties for maximizing information content in flight data is discussed, along with 
practical implementation issues. Results from flight tests of the X-43A hypersonic research 
vehicle (also called Hyper-X) are used to demonstrate the excellent modeling results 
obtained using this maneuver design approach. A detailed design procedure for generating 
the maneuvers is given to allow application to other flight test programs. 
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Nomenclature 

body-axis translational accelerometer measurements, ft/sec 2 
wing span, ft 

mean aerodynamic chord, ft 
mass moments of inertia, slug-ft 2 

frequency index 

number of inputs, or mass, slugs 

Mach number 

maximum 

minimum 

body-axis roll, pitch, and yaw rates, rad/sec 

dynamic pressure, lbf/ft 2 

root-mean-square 

wing reference area, ft 2 

excitation time length, sec 

airspeed, ft/sec 

coordinates of the center of gravity, in 
coordinates of the reference point, in 
angle of attack, deg 
sideslip angle, deg 

symmetric elevon, asymmetric elevon, and rudder deflections, deg 
phase angle, rad 
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I. Introduction 


O ne of the goals of current aerospace research is to develop vehicles that can maneuver at hypersonic speeds. To 
achieve that goal, it is necessary to conduct flight tests of hypersonic vehicles to determine their stability and 
control characteristics. Presently, it is not possible to sustain hypersonic flight at a given flight condition for any 
substantial length of time, because of limitations in hypersonic propulsion, thermal protection, and practical 
considerations such as airspace limitations and the size of available test ranges. Consequently, flight conditions are 
transient when flight testing hypersonic vehicles. In particular, the Mach number changes rapidly and significantly 
during hypersonic flight tests, which means that any maneuver must be executed efficiently in all axes over a very 
short time period. This is done so that the flight condition for each maneuver can be considered relatively constant, 
which simplifies the data analysis and modeling. If the maneuver takes too long, the Mach number changes 
significantly, and the vehicle aerodynamics, stability, and control change with it. 

Stability and control characteristics of a flight vehicle can be quantified by parameters appearing in the equations 
of motion, such as stability and control derivatives in the case of a linear dynamic model 12 . Values of these 
parameters can be found from flight data using various forms of least squares fitting of the model output (which 
depends on the values of the model parameters to be estimated) to measured data 1 . For hypersonic vehicles (and 
many experimental subsonic aircraft as well), flight test maneuvers are executed by applying designed input time 
series to the control surfaces, using an onboard flight computer 3 . These inputs are perturbations about the nominal 
or trim control surface deflections. 

The requirements for an effective maneuver to determine vehicle stability and control characteristics at 
hypersonic speeds are as follows: 

1 ) Multiple simultaneous inputs - For maximum flight test efficiency, and to ensure that the vehicle response 
in all axes is recorded at approximately a constant flight condition, all control surfaces on the vehicle should be 
moved simultaneously. This approach saves time compared to the conventional approach of moving one 
control surface at a time in sequence. 

2) Orthogonal inputs - With all inputs moving simultaneously for flight test efficiency, it is best if the inputs 
are mutually orthogonal. Mutual orthogonality means the inputs are the most unlike each other that they could 
possibly be. Since each input is doing something completely different, the effect of each individual control can 
be identified as accurately as possible. Considering that the data analysis could be done in either the time 
domain or the frequency domain, it would be best if this orthogonality could be enforced in both domains. 

3) Wide-band frequency content - The best inputs for exciting the vehicle dynamic response would have 
frequency content at or near the modal frequencies of the vehicle. However, these modal frequencies are only 
known approximately - if they were known exactly, there would be no need for flight test. A robust input 
design has wide-band frequency content to cover many possibilities. As long as the frequency content of each 
input is wide-band, the vehicle dynamic response will be excited, regardless of exactly where the modal 
frequencies lie. Another practical reason for wide-band inputs is that a single wide-band input design can be 
applied repeatedly at many different flight conditions, without sacrificing effectiveness. The practical result is 
a savings in onboard computer memory required, and a simplified input design. 

4) Small excursions - The inputs must excite the vehicle dynamic motion, but it is best if this motion can be 
limited to small excursions from the target flight condition, so that the vehicle stability and control can be 
adequately described using linear models. For hypersonic flight tests, the requirement for small excursions 
also arises from the desire to fly the vehicle on a specific planned trajectory. Even small deviations from a 
planned trajectory can cause large uncertainty in the landing or splashdown location. 

5) Perturbation inputs - To achieve small excursions about the target flight condition and keep the vehicle on 
its intended trajectory, the inputs for the maneuver must be perturbations about the nominal or trim control 
surface deflections. This means that each input must begin and end at zero amplitude, and oscillate around 
zero amplitude. 

In this paper, a method for designing inputs that meet the above requirements is described in detail. This 
method, first described in Ref. 3, has been applied successfully to the NASA X-43A (Hyper-X) hypersonic flight 
tests 4 ' 5 . The form of each input is a sum of sinusoids with unique frequencies and phase shifts. Multiple inputs are 
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designed to be mutually orthogonal in both the time domain and the frequency domain, and are optimized for 
maximum data information content in multiple axes over a short time period, while minimizing excursions from the 
target flight condition. The only a priori information required for the multiple-input design is an estimate of the 
approximate frequency band for the vehicle dynamics, and approximate relative control effectiveness for proper 
relative scaling of the input amplitudes. 

The next section describes the input design procedure in detail, along with some analytical results and practical 
considerations. This is followed by a description of the X-43A hypersonic research vehicle and flight test, and 
modeling results from the X-43A flight data. 

II. Input Design Method 


A. Multiple Input Design 

The general idea is to excite the vehicle using inputs with frequency content near where the modal frequencies of 
the vehicle dynamic response are expected to lie, using an input form that is efficient for parameter estimation. Each 
input is a sum of sinusoids with unique frequencies that cover a frequency band of interest, similar to frequency 
sweeps 1 , and the amplitudes of the sinusoidal components are chosen to achieve a specific power distribution. The 
wide-band frequency content of the inputs is important because there is naturally some uncertainty as to the exact 
modal frequencies for the vehicle in flight. Wide-band inputs provide robustness to this uncertainty. 

The following material describes in detail how to design multiple orthogonal phase-optimized multi-sines and 
why this particular input form is efficient for estimating stability and control parameters in dynamic models. 

Each perturbation input u - , applied to the j th individual control surface, is comprised of a set of summed 
harmonic sinusoids with individual phase shifts fa , 


Z , . i '2.71 k t 

A k sin | + < 


j = 1,2 


( 1 ) 


where M is the total number of available harmonically-related frequencies, T is the time length of the excitation, and 
A k is the amplitude for the k th sinusoidal component. The variable t represents a vector of A discrete time 

points, t = [?(0) / (I ) ... t (A -l)]^ , and Uj represents the vector of corresponding amplitudes for the / h input, 

Uj = | uj (0) Uj (l) ... Uj (A - 1) j . Each of the m inputs is comprised of selected components from the pool 
of M sinusoids with frequencies co k = Ink/T , k = l,2,...,M . 

If the phase angles </> k in Eq. (1) were chosen at random on the interval \-7t,jz\ rad, then in general, the various 
harmonic components would add together at some points to produce an input Uj with relatively large amplitude 

excursions. This is undesirable, because it can result in the dynamic system being moved too far from a reference 
condition selected for the experiment. To prevent this, the phase angles <j) k for each of the selected harmonic 
components are chosen to minimize relative peak factor RPF , defined by 13-6 
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Relative peak factor is a measure of the efficiency of an input for parameter estimation purposes, in terms of the 
amplitude range of the input divided by a measure of the input energy. The relative peak factor is scaled so that any 
individual sinusoidal component (such as any one of the summands in Eq. (1)) has RPF = 1 . Low relative peak 
factors are desirable for estimating dynamic model parameters, because the objective is to excite the dynamic system 
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with a variety of frequencies without driving it too far away from the nominal operating point. The target value for 
optimizing RPF of an input with many harmonic components (cf. Eq. (1)) is 1, to match the RPF for a single 
harmonic. In some cases, the RPF for a phase-optimized composite input can be slightly less than 1, as will be 
shown later. Inputs with low relative peak factors are efficient in the sense of providing good input energy to excite 
the dynamic system over a selected frequency band, with low amplitudes in the time domain. 

For a composite signal with more than one sinusoidal component, as in Eq. (1), the goal is to design the input for 
minimum RPF. This is done by adjusting the phase parameters </> k for the components of each input to minimize 
RPF. The resulting optimization problem is non-convex; however, a simplex algorithm 7 can be applied to find a 
solution. 

The integers k specifying the frequencies for the / th input u ; are selected to be unique to that input, but are not 

necessarily consecutive. A good approach for multiple inputs is to assign integers k to each input alternately. This 
is illustrated in Figure 1 for the X-43A flight test maneuver design. In this case there were 3 inputs, elevator, 
aileron, and rudder, and a total of 42 frequencies (M = 42) . The frequencies were interleaved among the three 

inputs to achieve wide-band frequency content for 
each input. This provided robustness to uncertainty of 
how each control would excite the dynamic modes of 
the vehicle. Furthermore, since each input has wide- 
band frequency content, a single input design can be 
applied at various flight conditions, which simplifies 
the flight test and requires a small amount of onboard 
computer memory. The power spectra for elevator and 
rudder were modified so that more excitation power 
was applied at middle frequencies where the dynamic 
modes excited by these controls were believed to be. 
The power spectra shown are normalized, so the 
effects of individual control surface amplitudes are 
excluded. This means that for each input, the sum of 
all the spectral line ordinates (sum of the heights of the 
bars for each color) is 1 . 

When the frequency indices k that are selected for 
each input Uj in Eq. (1) are distinct from those chosen 

for the other inputs, then the frequency content of each u ; consists of distinct spectral lines in the frequency 

domain, as can be easily seen in Figure 1. Therefore, the vectors of Fourier transforms for the inputs as a function of 
frequency have inner products equal to zero. In this sense, the inputs are mutually orthogonal in the frequency 
domain, because each input contains frequencies that no other input has. 

In the time domain, a sum of sinusoids is orthogonal to any other sum of sinusoids with harmonically-related 
frequencies, regardless of the constant phase shift of each sinusoidal component. For example, consider two inputs, 
each containing a single distinct harmonic sinusoidal component, 
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Figure 1 Multiple orthogonal phase-optimized 
multi-sine input spectra, X-43A 
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where k x and k 2 are both integers. Then using the discrete-time notation q = iAt and T = N At , where At is a 
constant time step, the inner product of these inputs is 
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Using the trigonometric identity 


sin (u^sin (v) = ^cos(n — v) — cos(u + v)^| 


(5) 


Eq. (4) is changed to 


T 1 V’ 

“l“2=-2.' 

A i = 0 


N-l 


2n(k\ -k 2 )i 


N 


N-l 


+ <h-<h.\ ~ - Yj cos 

1 i= o 


2 n(k\ +k 2 )i 


N 


+ ( k +< h. ~ 0 


(6) 


Each summation in Eq. (6) is zero because the sum of a harmonic sinusoid over the base period T is zero, 

^ (2 7tki \ 

- + </> = 0 , for any integer k and any constant phase angle <j) . This shows that harmonic sinusoids are 
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orthogonal in the time domain. For more than one sinusoidal component in each input, the analysis is similar. So, 
inputs assembled as in Eq. (1) are orthogonal in the time domain. 


An objective for the experiment design is to excite the aircraft dynamics in a short time period by moving 
multiple control surfaces simultaneously. Since more than one surface is being moved, it is advantageous for the 
modeling if the Uj vectors applied to each control surface 


are mutually orthogonal. This helps the parameter 
estimation by completely de-correlating the inputs, which 
improves the accuracy of control effectiveness estimates. 
Using the design method described here, it is possible to 
make all of the u j mutually orthogonal in both the time 

and frequency domains, using inputs designed for 
minimum relative peak factor. This gives the analyst the 
flexibility to use either time domain or frequency domain 
parameter estimation methods while retaining the desirable 
feature of mutually orthogonal inputs. 

Figure 2 shows perturbation inputs applied at Mach 6 
on the X-43A. These inputs are mutually orthogonal in 
both the time and frequency domains. The inputs were 
computed from Eq. (1), using the spectra shown in 
Figure 1, with phase angles fa optimized for minimum 
relative peak factor. Table 1 contains all the information 
for the input design at Mach 6. Because of the various 
frequencies and phase angles, and the small amplitudes of 
the input perturbations, applying these inputs 
simultaneously to the vehicle produces a dynamic response 
similar to what might be seen in flight through light to 
moderate turbulence. Consequently, the vehicle stays on 
its nominal trajectory, but jiggles around that. 



Figure 2 Multiple orthogonal phase-optimized 
multi-sine inputs, X-43A, Mach 6 


For each input, the power spectrum can be tailored by selecting the A k to distribute power over the spectral 
components. To achieve a uniform power distribution. 


A k 
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where n is the number of sinusoidal components included in the summation of Eq. (1), and A is the amplitude of 
the composite input Uj . Therefore, with uniform power distribution, selection of the A k reduces to selecting a 

single input amplitude A . 

To summarize, the multiple input design procedure is as follows: 

1) 1. Select the time period T for the excitation, which determines the smallest harmonic frequency resolution 
Af = \/T Hz and the limit on the minimum frequency, f min > 2/T Hz . Often this is dictated by practical 
considerations such as length of time available at a particular flight condition. 

2) Select the frequency band of the dynamic system for the excitation frequencies, \ f m j n ,f max ] Hz. This 
corresponds to the frequency band where the dynamic response of the system is expected to lie. The 
frequencies are equally spaced by A f on the interval [/„„■„ , f max ] ■ The total number of frequencies is 

M - fix\(f max - fmin ) /H / ] + 1 , where fix indicates rounding to the nearest integer toward zero. 

3) Assign approximately an equal number of indices k from the set { 1 , 2, . . . , A7| to each input by alternating 

each consecutive frequency among the multiple inputs. This approach produces lower relative peak factors for 
the individual inputs, and also ensures that each input has frequency content distributed evenly across the 
frequency band [/ m /„ , f m ax ] • Different assignments of the frequency indices could be made for other reasons. 

For example, lower frequency indices might be assigned to an input that is known to effectively excite a low 
frequency mode, or particular frequency indices might be omitted to avoid exciting an undesirable structural 
response or to simplify the input design. Each frequency index can be assigned to only one of the inputs, to 
preserve mutual orthogonality of the inputs in both the time and frequency domains. 

4) Generate the input u ; for each of the m controls, using Eq. (1) and assigning the starting values for the phase 
angles <f> k at random from a uniform distribution on the interval \-n, 7i\ rad. 

5) Choose the input power spectrum by specifying the amplitude A k for each harmonic component, or use 
Eq. (7) to implement a uniform power distribution. 

6) Apply a simplex optimization algorithm 8 to adjust the (j) k for each u f to achieve minimum relative peak 
factor for that input. The simplex optimization technique does not require gradients. 

7) For each input, do a one-dimensional search to find a constant time offset for the components of each input 
Uj , so that the input begins and ends at zero amplitude. This is equivalent to sliding the input along the time 

axis until a zero crossing is placed at the origin of the time axis. The appropriate phase shift is added to each 
component phase angle (f> k . Note that to implement a constant time shift to all the components, the phase 
offset for each component will be different, because each component has a different frequency. Since the 
components of each Uj are harmonics of the base frequency with period T , if all the component phase angles 

(f> k are shifted so that the initial value of the input is zero, then the final value of the input at time T will also 
be zero. The power spectrum, input orthogonality, and relative peak factor are all unaffected. 

8) Return to step 6 unless the relative peak factor has reached a pre-defined goal value or a maximum number of 
iterations has been reached. For practical work, the relative peak factor goal is set at 1.01 and the maximum 
number of iterations is set to 50. 
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B. Practical Issues 

All of the steps detailed above are implemented in MATLAB® code that is part of a software package called 
SIDPAC (System IDentification Programs for AirCraft). SIDPAC is bundled with Ref. 1. SIDPAC also contains 
many other useful programs for experiment design, data analysis, modeling, and simulation. The program that 
implements the above steps is called mkmsswp.m. 

Sometimes repeating the input design with 
different random starting values for the fa (cf. step 4 
above) can lead to an improved (lower) optimized 
relative peak factor solution. This behavior is 
common for high-dimensional non-convex 
optimization problems such as this one. In practice, 
this just means re-running the input design program, 
because each run uses a different realization of the 
random initial values for fa . 

To gain an appreciation of the non-convexity of 
the optimization problem, Figure 3 shows the relative 
peak factor as a function of phase angles for two 
components in the summation of Eq. (1), using 7= 1 5 
sec, unit amplitudes, and frequency indices k=2 and 
k= 4. It is clear from the figure that the optimization 
problem is not convex, and also that there are several 
phase angle solutions which are equally good, or 
nearly so. The data used to make Figure 3 can be used 
for a global exhaustive search for the minimum relative peak factor in this simple case. To generate this data, the 
relative peak factors were computed for values of the phase angles on a 2-dimensional grid with intervals of 
0.0175 rad (1 deg), over the range [ 0, n ] rad ( [ 0, 1 80 ] deg ) . The minimum relative peak factor of all these 

computed values was found to be 1.102. Using the technique described above to optimize the phase angles, the 
relative peak factor achieved was 1.106. This demonstrates that the technique developed here found an input with 
relative peak factor very close to the global optimum, for this simple case. The same exercise was then repeated for 
three harmonic components, with frequency indices k=2, 4, and 6, and the same phase angle intervals, in three 
dimensions. In this case, the global exhaustive search gave a minimum relative peak factor equal to 1.002, and the 
optimization technique described above produced an input with relative peak factor equal to 1.002. These 
investigations suggest that the optimization procedure does an excellent job of designing inputs with relative peak 
factor very close to the global optimum. 

Step 7 of the input design procedure ensures that each designed input will be a perturbation, so that the designed 
input can be added to whatever value the control may have for another reason, e.g., trim or maneuvering. The 
iteration in this step helps to reduce the relative peak factors by adding perturbations to the phase angles <f> k (from 
the phase changes added to make the input start and end at zero), which helps the optimization algorithm solve the 
non-convex optimization problem. 

It was found empirically that if the set of selected harmonic frequency indices k selected for a particular input 
consisted of an integer greater than 1, along with 2 or 3 multiples of that integer (e.g., k= 2, 4, 6, or k= 5, 10, 15, 20), 
the phase angles could be optimized (in step 6) so that the relative peak factor for that input was very close to 1, and 
sometimes less than 1 . 

Time length of any input can be extended or reduced, without affecting the orthogonality properties or the 
relative peak factors, by simply re-computing each input from Eq. (1) using the same optimized phase angles fa 
and amplitudes A k , but a different value of T . This is equivalent to stretching or shrinking each harmonic 
component. Of course, if T is changed, the frequency content of the input changes. As T increases, the 
frequencies in the input shift lower, with decreased frequency spacing. This is demonstrated in Figure 4 for a single 
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optimized input form with frequency indices k = 2, 4, 6, 8 . The phase angles of the harmonic components for each 
of the inputs shown in Figure 4 are the same, and the relative peak factor of each input is 0.995. 
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Figure 4 Effect of changing excitation time T 

It is possible to choose non-integer values of k when assembling the inputs using Eq. f 1). This allows arbitrary 
spacing of the frequencies for each input. As long as unique frequencies are assigned to each input u .■ , mutual 

orthogonality in the frequency domain remains intact. However, Eqs. (3)-(6) and the associated discussion show 
that inputs with non-integer values of k will not be mutually orthogonal in the time domain unless the sum and 
difference of the frequency indices are all equal to an integer. It may still be advantageous to use non-integer values 
of k in some cases, because even though the inputs will not be orthogonal in the time domain, their correlation will 
still be low, and this situation is still quite acceptable for good modeling results. Furthermore, if the vehicle has a 
feedback control system operating during the experiment, then the time-domain orthogonality of the inputs is mined 
anyway, by the actions of the feedback control. Assuming the resulting input correlation is not too high (absolute 
value of the correlation must be below 0.9), the modeling results will not be adversely affected. 

Although there are methods for optimizing the frequency spectrum of inputs for parameter estimation (i.e., the 
choice of A k in Eq. (1)), all of them require some knowledge of the vehicle dynamics, usually in terms of a nominal 
dynamic model with parameter values. For hypersonic flight tests, there is questionable value to such methods, 
because there is typically significant uncertainty associated with pre-flight predictions of hypersonic vehicle 
dynamics. Under these circumstances, it is best to define the frequency spectra for all inputs to be relatively flat 
across the selected frequency band, so that the aircraft dynamics will be sufficiently excited, regardless of where the 
modal frequencies happen to be located within the frequency band. 

The next two sections describe how the input design method explained above was applied to flight testing the 
X-43A (Hyper-X) hypersonic vehicle. 
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III. Research Vehicle and Experimental Data 


100,000 



40,000 


- Hyper-X free flight 


Descent 


On November 16, 2004, flight 3 of the NASA X-43A hypersonic research vehicle, or Hyper-X, was successfully 
conducted over the Pacific Ocean off the California coast. The flight included a successful test of the scramjet 

propulsion system at approximately 
Mach 10. Following the engine test, the 
engine cowl door was closed, and a 
series of maneuvers were executed for 
aerodynamic parameter estimation as the 
vehicle descended and Mach number 
decreased. Maneuvers were executed for 
flight conditions at Mach 8, 7, 6, 5, 4, 
and 3. The vehicle flew as a glider 
(engine off) throughout these maneuvers. 
Figure 5 shows the X-43A mission 
profile. Figure 6 is an artist’s rendering 
of the X-43A research vehicle in flight. 

A 3 -view drawing is shown in Figure 7. 


■ Scramjet 
engine start 

Hyper-X booster 
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Air launch 


Booster burn-out 


(Over water) 


Experiment 1 
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Distance 


Figure 5 Hyper-X mission profile 



Figure 6 X-43A research vehicle 



Figure 7 X-43A research vehicle 3-view 


At the start of the mission, the X-43A research vehicle was mounted to the nose of a modified Pegasus booster. 
This configuration, called the Hyper-X stack, was carried aloft by the NASA B-52B aircraft in a manner similar to 
the X-l flight tests. Figure 8 is a picture of the Hyper-X stack being carried to altitude for flight 3. At 40,000 ft, the 
Hyper-X stack was dropped, and the Pegasus booster fired about 5 seconds later. The Pegasus booster carried the 
X-43A to the target flight condition for the scramjet engine test, which was approximately Mach 10 and 100,000 ft 
altitude. At that point, the X-43A separated from the Pegasus 
booster, and began its free flight. After stabilization at the target 
flight condition, the engine cowl was opened and the scramjet 
engine test was conducted. Then, the engine cowl door closed, 
and the X-43A executed an energy reduction maneuver to arrive 
at the Mach 8 flight condition. Perturbation inputs for 
aerodynamic parameter estimation, control law validation, and 
aerodynamic prediction were then run automatically, and these 
were repeated at Mach 7, 6, 5, 4, 3, and 2 during the vehicle 
descent. The form of these inputs (shown in Figure 2) was the 
same at each Mach number, but the amplitude of the inputs 
changed with flight condition. The vehicle then splashed into the 
Pacific ocean and was not recovered. Note that the data at Mach 
2 was not usable for parameter estimation, because at that flight 

condition, the vehicle encountered wind shears and gusts which Figure 8 Hyper-X stack being carried to 

were unmeasured and therefore corrupted the parameter altitude by the NASA B-52B 
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estimation results. Consequently, flight 3 results did not include the Mach 2 flight condition. 

Ref. 8 includes complete information on the geometry and mass properties of the X-43A research vehicle. Table 
2 gives a summary of the constant mass and geometry values used for analysis of the descent data. The mass 
properties reflect the fully loaded X-43A research vehicle minus the mass of fuel and coolant used during the 
scramjet engine test. 

Flight data came from the Best Estimated Trajectory (BET), documented in Refs. 8 and 9. The BET data set used 
the best data from several recorded sources, corrected for known time skews and estimated winds aloft. The BET 
data were also corrected for systematic instrumentation errors using a data compatibility analysis. Sample rate for 
the BET flight data was 50 Hz, corresponding to a sampling interval of 0.02 sec. 

All parameter estimation was done for the aerodynamic forces and moments acting at the model reference point, 
which is the reference point for the aerodynamic database. This was done to facilitate comparisons between stability 
and control derivatives from the aerodynamic database and estimates based on flight data. The same approach was 
used for the flight 2 results documented in Ref. 4. Details of the location of the model reference point are available 
in Ref. 8. Coordinates of the reference point are given in Table 2, for a standard right-handed coordinate frame using 
fuselage station, butt line, and water line coordinates, as defined in Ref. 8. The aerodynamic reference point is about 
5.5 feet aft of the nose, on the vehicle centerline. 

IV. Flight Test Results 

The descent maneuvers were designed using the procedure described in Section II. Multiple orthogonal 
phase-optimized multi-sine inputs were applied to all three control surface inputs - symmetric elevon (elevator), 
asymmetric elevon (aileron), and rudder. Since flight conditions change rapidly during the descent, the inputs were 
designed to be applied simultaneously, so that a full set of non-dimensional longitudinal and lateral stability and 
control derivatives could be estimated at each flight condition. The flight testing also included multi-step 2-1-1 
waveforms applied to each of the control surface inputs in sequence, following the multi-sine inputs. Flight data 
from the 2-1-1 maneuvers were used to test the prediction capability of the models identified from the multi-sine 
inputs. Accurate prediction for maneuvers with dissimilar inputs is a strong indicator of good modeling results. 

Figure 2 shows the multiple orthogonal phase-optimized multi-sine perturbation inputs that were applied at 
Mach 6. The same input forms were used at every flight condition from Mach 8 to Mach 3, but the amplitudes of 
the inputs were adjusted to account for the changing dynamic pressure and control surface effectiveness. This 
approach was taken because of a practical limitation on flight computer memory. The wide-band frequency content 
of the inputs made it possible to satisfy this practical constraint without compromising the effectiveness of the 
excitation. 

As noted earlier, orthogonal inputs are advantageous for parameter estimation experiments, because the 
aerodynamic dependencies can be assigned more accurately when all inputs are doing something different to excite 
the dynamic response of the vehicle. The orthogonality of the designed inputs was disrupted during the flight testing 
of the X-43A, due to the action of the feedback control, which could not be turned off. However, starting out with 
orthogonal inputs prevented the feedback control from correlating the inputs to the point where the modeling results 
would be adversely impacted. In effect, since the inputs were orthogonal (i.e., perfectly decorrelated) to start with, 
the degradation by the feedback control only reduced them to very good (i.e., very low correlation). 

Amplitudes for the three inputs were selected based on test runs of a nonlinear simulation of the Hyper-X, which 
used an aerodynamic database 10 built from extensive wind tunnel tests and CFD, and included the feedback control 
system. The aim was to scale the inputs to achieve response amplitudes of approximately +3 deg in angle of attack 
and sideslip angle perturbations for each maneuver, based on the simulation. Input amplitude adjustments were 
done by trial-and-error, because the feedback control system distorted the perturbation inputs for parameter 
estimation in a manner that depended on the vehicle response. Note that the form of the perturbation input 
commands did not change - this form was optimized by the procedure described earlier. The trial-and-error was 
only for sizing the amplitudes of the three inputs at the various flight conditions from Mach 8 to Mach 3. 

Standard routines from the software toolbox called System IDentification Programs for AirCraft ( SIDPAC )’ 
were used to identify the appropriate model structure for the aerodynamic forces and moments, and to estimate the 
parameters in the identified model structures. Details on how this was done can be found in Refs. 1, 3, and 4. 
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For all flight conditions from Mach 8 to Mach 3 during the descent, a linear aerodynamic model was found to be 
adequate to characterize the measured data. The aerodynamic model equations were 


Ct — C j cc + C i S p + C j 

L L a L S e e L o 

(8) 

Cd = Cd u a + Cog e 3 e + Q) 0 

(9) 

(- 'm ~ C m a ^ C m q ^ ^m 0 

(10) 

Cy = C Yf , P + C Y g r <5 r + C Yo 

(ID 

C l =C lpP + C h a 5 a + C l Sr S r + C l a 

(12) 

c„ = c llfj p + c nSa s a + c ngr S r + C„ o 

(13) 


where all state and control variables are perturbations from a reference condition, defined at the beginning of each 
maneuver. 

Lift and drag coefficients were used instead of the body-axis X and Z components of aerodynamic force, so that 
the flight results could be compared directly to values from the pre-flight aerodynamic database 10 . Aerodynamic 
moment coefficients were modeled at the aerodynamic reference point used for the aerodynamic database, again to 
facilitate comparisons. 

For all flight conditions studied, the pitch rate damping was the only rate derivative term that had a significant 
effect on the X-43A aerodynamics. Model terms associated with other rate derivatives were dropped during model 
structure determination, because the models given in the above equations almost completely characterized the 
variations in the associated force or moment coefficient. This simply means that rate damping terms were not 
necessary to model the lift, drag, side force, or lateral moment coefficients, so the damping terms were omitted from 
those models. Equivalently, the damping parameters were zero for these force and moment coefficients, based on 
the flight data. 

The aerodynamic model parameters in Eqs. (8)-(13) can be estimated from flight data using linear regression or 
output-error maximum likelihood in the time-domain or frequency-domain, among other approaches 1 . Software that 
implements each of these methods is included in SIDPAC, and these tools were applied to the flight data. The most 
consistent results, which also exhibited the best prediction capability, came from equation-error in the frequency 
domain. Consequently, all of the parameter estimation results shown in this work were computed using 
equation-error in the frequency domain. Detailed information on this approach can be found in Refs. 1, 3, and 4. 

The equation-error method calculates aerodynamic parameter estimates that minimize squared errors between 
values of the non-dimensional force and moment coefficients determined from measured flight data, and model 
values computed from Eqs. (8)-(13). For a glider (e.g., the X-43A in descent), non-dimensional force and moment 
coefficients are determined from measured flight data using the following equations: 


_ _ma 
c x = ~ l a ~ 


ij S 


Cy = 


Y ~ — 


qS 


(14) 


( 15 ) 
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( 16 ) 


mci- 

-Z =~'~N = ~ 


qS 


C/ = —C z cos a + C x sin a 
C D = -C x cos a -C z sin a 


C,= 


i -~ 


qSb 


I xz( ( Jz ly) 

p -j Hpq +r ) + j — V 


c = y ~ 

c-m — c — 

qSc 


q +^Y~^~P r + -j 1 -(p 2 _r2 ) 


C n =-^~ 

qSb 


r-^(p-qr)+^—^-pq 

*z *z 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


Substituting measured data into the right sides of Eqs. (14)-(21) results in (V values of the non-dimensional force 
and moment coefficients, where N is the number of data points for the maneuver. Using these values and measured 
states and controls in Eqs. ( 8)-( 1 3) results in an over-determined set of equations for the unknown aerodynamic 
parameters, which can be solved with standard least-squares methods 1-3 ' 4 ' 7 . 

Estimates for stability and control derivatives associated with the lift coefficient are plotted in Figure 9 as a 
function of Mach number. The triangles represent flight estimates, and the circles are values computed from table 
look-ups in the pre-flight aerodynamic database. Error bars on all estimates indicate +2 standard errors, 
corresponding to 95% confidence, assuming a Gaussian distribution. Error bars for the aerodynamic database values 
represent an estimated uncertainty, computed based on the uncertainty model documented in Ref. 1 1 . Error bars for 
the flight estimates were statistical values based on the flight data. 

Numbers on vertical axes for the parameter estimate plots have been removed, because the numerical results are 
restricted to U.S. government employees and U.S. government contractors only. However, the scales have not been 
altered, so the relative positions of the estimates and sizes of the error bars are real. 

Similar plots for pitching moment coefficient parameters and drag coefficient parameters are shown in 
Figures 10 and 11, using the same presentation format. Side force, yawing moment, and rolling moment coefficient 
parameters appear in Figures 12, 13, and 14, respectively. Inspection of these plots shows that the flight estimates 
were generally in good agreement with values obtained from the pre-flight aerodynamic database. This is indicated 
by flight estimates usually lying within the uncertainty bounds of the corresponding value from the aerodynamic 
database. The flight estimates had small error bars, indicating accurate parameter estimates, which is the direct 
result of effective experiment design, including both instrumentation and the input design for the maneuvers. The 
standard errors for the flight estimates were less than 5 percent of the estimated values, for all parameters except the 
pitch rate derivative C m . No adjustments were made to the flight estimates or the flight data analysis procedure to 

improve the match to the aerodynamic database. Flight estimates and their error bounds were computed first, and 
separately, then simply plotted with the values from the pre-flight aerodynamic database as Figures 9-14. 

The pitch damping C m in Figure 10 shows what looks like a significant mismatch between flight estimates and 
the aerodynamic database. However, the C m parameter multiplies the quantity qc /2V in the model equations, 
cf. Eq. (10). At Mach 5, for example, the value of qc /TV is 0.0006 for q = 0.5 rad/sec, which means that the value 
of C m must be quite large in order for the damping term C m to be a non-trivial contributor in the model of 
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Eq. (10). In fact, the estimates of C m shown in Figure 10 were small in magnitude, so that the mismatch seen for 
the C m estimates makes very little difference in the pitching moment computed from Eq. (10). The possibility of a 

numerical problem in the parameter estimation arising from the relatively small value of qc /2V compared to the 
other explanatory variables in the model was examined and ruled out. 

Modeling results from flights 2 and 3 could not be compared for maneuvers executed at similar Mach numbers, 
because the flight conditions were different. Consequently, it was more meaningful to compare flight 3 estimates to 
values extracted from the aerodynamic database at the same flight conditions, rather than to compare flight estimates 
for flights 2 and 3. The maneuvers at similar Mach numbers on flights 2 and 3 had dissimilar values of important 
flight condition parameters such as nominal angle of attack and elevator deflection. Agreement of flight 2 parameter 
estimation results with values from the pre-flight aerodynamic database (cf. Ref. 4) was very similar to that shown 
in Figures 9-14 for flight 3. 

Figure 15 shows a typical model fit to the flight data. This particular example is the model fit to the pitching 
moment coefficient at Mach 5 on flight 3. The excellent fit quality is typical of the other cases as well. The upper 
plot in Figure 15 shows the model fit to the data in the frequency domain, the middle plot shows the model fit to the 
data in the time domain, and the lower plot shows the difference between the two traces in the middle plot, which is 
the residual in the time domain. The residual shows no deterministic content, indicating that the model has captured 
all of the deterministic variation. Note that the time-domain model includes a time-domain least-squares estimate of 
bias and linear trend with time in the measured output, because these terms are omitted when modeling in the 
frequency domain. 

In Figures 16 and 17, the model identified using flight data from a maneuver designed as described above was 
used to predict the pitching moment coefficient at a similar Mach number (approximately Mach 5), using a different 
input, namely a 2-1-1 multi-step. Figure 16 shows the Mach number variation over the time corresponding to the 
parameter estimation input (done first), followed by a control system validation sweep at higher frequencies, then 
finally the 2-1-1 prediction maneuver. The upper plot in Figure 17 shows an expanded plot of the 2-1-1 command to 
the symmetric elevon appearing at the end of the sequence in Figure 16, along with the measured symmetric elevon 
position. The commanded symmetric elevon was modified by feedback control. The middle plot in Figure 17 
shows the time-domain prediction of pitching moment coefficient and the lower plot is the residual, or the difference 
between the two traces in the middle plot. The quality of the model fit to the prediction data is similar to that seen in 
Figure 15, indicating that the identified model is a good predictor for a dissimilar maneuver at similar flight 
conditions. Good prediction capability is a very strong indicator of an accurate model. The quality of the prediction 
result shown in Figure 17 was typical of the prediction results for other flight conditions. 

Slight mismatches toward the end of the prediction maneuver in Figure 17 might be attributable to changes in the 
vehicle aerodynamics due to rapidly changing Mach number during the descent. Figure 16 shows that the Mach 
number variation over the period of time from the parameter estimation input to the 2-1-1 prediction maneuver is 
approximately -0.75. The model identified from data near the start of the time period may be slightly inapplicable 
for the prediction maneuver executed near the end of the time period. 

V. Concluding Remarks 

A multiple-input design technique for characterizing the stability and control of hypersonic vehicles was 
described and demonstrated. The technique can be used to design multiple inputs that are mutually orthogonal in 
both the time and frequency domains. Mutual orthogonality of the inputs is important for accurate modeling results, 
and also improves flight test efficiency and practicality, because orthogonal inputs can be applied simultaneously. 
This allows all relevant data to be collected at an approximately constant flight condition, simplifying the data 
analysis and modeling. Since the maneuvers are designed for selected frequency ranges, there is no need for a priori 
models in the design process, and the same input forms can be used over a wide range of flight conditions. In 
addition, the inputs are optimized for minimum amplitude excursions in the time domain for a given input energy. 
Power spectra for the inputs can be chosen to target or exclude specific frequencies or frequency bands. These 
features make the inputs practical and useful for hypersonic flight tests, where the excitation must have short 
duration and low output response amplitudes, and modal frequencies vary with rapidly changing flight conditions. 
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The practical application of the multiple input design technique was demonstrated by estimating stability and 
control parameters for the X-43A hypersonic research vehicle from flight data. Maneuvers were flown using 
multiple orthogonal phase-optimized multi-sine inputs applied as simultaneous control surface perturbations at 
Mach 8, 7, 6, 5, 4, and 3 during the vehicle descent on flight 3. All model parameters, consisting of 
non-dimensional longitudinal and lateral stability and control derivatives, were estimated based on flight data from a 
single multi-axis maneuver executed at each Mach number. Multi-step inputs at nearly the same flight conditions 
were also flown to assess the prediction capability of the identified models. Parameter estimates were plotted as a 
function of Mach number, and compared with estimates from a pre-flight aerodynamic database, which was built 
from extensive wind-tunnel tests and computational fluid dynamics. Accuracy of the flight-determined modeling 
results was significantly better than the estimated uncertainty in the pre-flight aerodynamic database. 

Multiple orthogonal phase-optimized multi-sine inputs were shown to be good practical inputs for determining 
stability and control characteristics of hypersonic vehicles. These wide-band inputs provided sufficient excitation to 
the X-43A in spite of a feedback control system that distorted the desired input forms. The resulting flight data had 
sufficient information content for the extraction of highly accurate estimates of aerodynamic stability and control 
derivatives. Standard errors for the aerodynamic parameters estimated from flight data were less than 5 percent for 
all parameters except the pitch damping derivative. Flight data from multi-step 2-1-1 maneuvers at similar flight 
conditions were used to show that the models identified using data from the multiple orthogonal phase-optimized 
multi-sines had good prediction capability for maneuvers with dissimilar input forms. All of these results support 
high confidence in the models identified from flight data. 

In addition to assessing vehicle stability and control, aerodynamic parameters estimated from hypersonic flight 
data are useful for comparisons with results from wind tunnel tests and computational fluid dynamics (CFD). This 
helps to validate and refine aerodynamic prediction methods for designing future hypersonic vehicles. 
Flight-determined aerodynamic parameters can also be used to update aerodynamic databases for hypersonic 
vehicles, for improved simulation, control system design, and dynamic analysis. Of course, the multiple orthogonal 
phase-optimized multi-sine inputs can be applied very effectively to other dynamic systems as well, for many of the 
same reasons outlined here for hypersonic vehicles. 

Because of the success of the X-43A flight tests and the practical considerations discussed earlier, the input 
design method described in this paper is being applied to the upcoming hypersonic flight testing for stability and 
control characterization on the U.S. Air Force / DARPA FALCON project 12 . 
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Table 1 Multiple input design for the X-43A 
flight 3 descent, Mach 6, T =27 sec 


Input 

A (deg) 

A k (deg) 

k 

fa ( rad > 

RPF 

s e 

1.0 

0.3062 

4 

0.1356 

1.179 

0.3536 

7 

2.7368 

0.3536 

10 

3.5216 

0.3536 

13 

3.4517 

0.3536 

16 

1.2786 

0.3536 

19 

3.1393 

0.3536 

22 

4.3436 

0.3536 

25 

4.6801 

0.3536 

28 

2.7439 

0.3536 

31 

-0.1261 

0.3307 

34 

2.6789 

0.3062 

37 

2.4822 

0.2795 

40 

0.8446 

0.2500 

43 

4.3635 

S a 

0.4 

0.2673 

5 

3.1024 

1.215 

0.2673 

8 

3.1671 

0.2673 

11 

0.7765 

0.2673 

14 

4.1085 

0.2673 

17 

-1.0068 

0.2673 

20 

0.9306 

0.2673 

23 

-0.0231 

0.2673 

26 

4.4717 

0.2673 

29 

0.4577 

0.2673 

32 

0.8148 

0.2673 

35 

1.8706 

0.2673 

38 

0.2096 

0.2673 

41 

4.5988 

0.2673 

44 

-0.1935 



1.4697 

3 

3.8424 


5 r 

5.0 

1.6971 

6 

4.0608 

1.029 



1.6971 

9 

3.4827 




1.6971 

12 

0.7669 




1.6971 

15 

3.2681 




1.6971 

18 

2.8489 




1.6971 

21 

2.0128 




1.6971 

24 

0.6371 




1.6971 

27 

2.7851 




1.6971 

30 

-1.4702 




1.5875 

33 

3.6098 




1.4697 

36 

2.5226 




1.3416 

39 

-1.3918 




1.2000 

42 

1.2756 
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Table 2 Geometry and mass properties 
of the X-43A research vehicle for flight 3 descent 


c , ft 

12.00 

b, ft 

5.19 

S , ft 2 

36.144 

X ref , in 

66.240 

y ref , in 

0.000 

Z ref , in 

0.000 

x c % , in 

69.193 

y CK * in 

0.019 

Z CK > in 

-0.547 

mass m , slugs 

87.75 

I x , slugs-ft 2 

54.46 

l y , SlUgS-ft 2 

851.40 

I z , slugs-ft 2 

880.37 

I xz , slugs-ft 2 

23.51 
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Figure 9 Lift parameters 



Figure 11 Drag parameters 
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Figure 10 Pitching moment parameters 
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Figure 12 Side force parameters 
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Figure 13 Yawing moment parameters 
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Figure 14 Rolling moment parameters 
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Figure 15 Model fit to pitching moment data, 
Mach 5 
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Figure 1 7 Model prediction of pitching 
moment data, Mach 5 



Figure 16 Mach number variation, 
Mach 5 maneuvers 
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